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Abstract 

The  aim  of  this  work  is  to  develop  a  new  rigorous  model  to  study  the  structure  performance  relationship  of  solid  oxide  fuel  cell  (SOFC)  electrodes. 
A  new  two-dimensional,  geometrical  model  which  captures  the  inhomogeneous  nature  of  the  location  of  electrochemical  reactions  based  on  random 
packing  of  electronic  and  ionic  conducting  particles  has  been  developed.  The  results  show  that  the  concentration  of  oxygen  inside  the  cathode  in 
the  2D  model  is  not  only  a  function  of  the  electrode  depth  but  also  changes  along  the  width  of  the  electrode.  Furthermore  the  effect  of  composition 
of  the  electrode  on  the  length  of  active  three-phase  boundary  (TBP)  and  total  polarization  resistance  has  been  demonstrated.  A  parametric  study 
of  the  effect  of  the  conductivity  of  ionic  conductor  and  diffusion  coefficient  on  the  performance  of  the  electrode  has  been  given. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  introduction  of  composite  electrodes  for  solid  oxide 
fuel  cells  (SOFC’s)  has  opened  a  new  promising  perspec¬ 
tive  toward  development  and  commercialization  of  SOFC’s 
[1].  Despite  simple  single  phase  electrodes,  cermet  compos¬ 
ite  electrodes  are  mixture  of  electronic  and  ionic  conductors 
and  hence  the  electrochemically  active  region  is  not  limited 
to  the  electrode-electrolyte  interface  and  the  reaction  extends 
inside  the  electrode.  As  a  result  of  increase  in  electrochemically 
active  area,  cermet  composite  electrodes  improve  the  solid  oxide 
fuel  cell  performance  and  allow  larger  power  densities  to  be 
obtained  [2,3]. 

In  the  present  work  we  focus  on  modelling  the  electrochem¬ 
ical  processes  at  the  cathode  in  order  to  develop  a  better  under¬ 
standing  of  structure  performance  relationship  of  electrodes. 
Researches  have  been  conducted  in  the  past  to  study  the  effect 
of  material  and  preparation  methods  on  microstructure  and  per¬ 
formance  of  SOFC  electrodes  [4-7].  Although  a  wide  variety  of 
materials  have  been  proposed  and  examined,  lanthanum  stron¬ 
tium  magnite  or  LSM  (Lai_xSrxMn03)  is  considered  to  be  the 
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most  promising  electron  conductor  along  with  Yttria- stabilized 
zirconia  (YSZ)  as  an  ion  conductor.  Apart  from  relatively  high 
electronic  conductivity  and  low  ionic  conductivity,  LSM  has 
good  electrocatalytic  properties  for  the  reduction  process  that 
takes  place  in  the  cathode  and  its  thermal  expansion  coefficient 
is  very  close  to  that  of  YSZ,  which  makes  this  pair  a  popular 
combination  for  use  in  composite  electrodes  [8-10]. 

Many  parameters  influence  the  performance  of  SOFC  com¬ 
posite  cathodes.  Among  those  are  size  of  electronic  and  ionic 
conductor  particles,  composition  of  the  cathode,  porosity  and 
electrode  thickness.  Although  some  experimental  work  has 
been  done  to  determine  the  effect  of  those  parameters  [3,7,11], 
broader  study  of  SOFC  cathodes  through  detailed  simulation  to 
obtain  an  optimized  structure  seems  to  be  necessary. 

Few  models  have  been  proposed  by  the  previous  researchers 
to  describe  the  SOFC  cathode  operation.  Despite  minor  differ¬ 
ences  in  interpretation  of  reaction  kinetics,  most  of  these  models 
rely  on  the  homogeneity  of  composite  electrode  media  for  the 
electrochemical  reaction  and  transport  of  ion,  electron  and  mate¬ 
rial  [12-14]. 

Costamagna  et  al.  [12]  considered  whole  contact  area 
between  electronic  and  ionic  conductors  to  be  active  for  the 
electrochemical  reaction  and  calculated  the  active  area  using  per¬ 
colation  theory  and  particle  coordination  number.  They  showed 
that  simultaneous  electronic  and  ionic  conductivity  inside  the 
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Nomenclature 

C  constant  equal  to  1.883  x  10-2 

Co  concentration  of  oxygen 

D  diffusion  coefficient 

eff  effective  properties 

el  electronic  conductor 

/  constant  equal  to  F/RT  where  F  is  Faraday’s  con¬ 

stant  and  R  is  gas  constant 

i  current  density  generated  on  the  boundary 

(Am-2) 

z'o  exchange  current  density  (A  m-2) 

io  ionic  conductor 

J  current  density  inside  ionic  and  electronic  media 

(Am-2) 

M  molecular  mass  of  gaseous  species  (gmol-1) 

P  pressure  (Pa) 

T  temperature  (K) 

Greek  letters 

a  charge  transfer  coefficient  in  Butler- Volmer 

equation  (=0.5) 

s  porosity  of  the  cathode 

< p  potential  (V) 

0  overpotential  (V) 

p  Resistivity  (Q  m) 

a  characteristic  length  of  atoms  in  Eq.  (9) 

r  tortuosity  of  the  cathode 

£?d  diffusion  collision  integral  in  Eq.  (9) 


electrode  is  not  likely  to  happen  for  particles  of  the  same  size 
when  concentration  of  one  of  the  species  is  less  than  about 
30  vol.%  and  depending  on  the  relative  size  of  ionic  and  elec¬ 
tronic  particles  they  obtained  the  minimum  polarization  resis¬ 
tances  for  the  compositions  of  the  electrode  around  40  vol.%. 
Chen  et  al.  [15]  proposed  a  new  kinetic  model  for  electrochem¬ 
ical  reaction  and  showed  that  polarization  resistance  assumes 
a  relatively  constant  value  over  the  wide  composition  range  of 
about  37-63  vol.%  of  YSZ.  They  showed  that  the  electronic 
current  generation  is  much  higher  near  the  electrode-electrolyte 
interface  and  attributed  this  observation  to  the  low  ionic  conduc¬ 
tivity  of  YSZ  [13]. 

Sunde  [16]  on  the  other  hand  used  Monte  Carlo  method  to 
show  that  while  changing  the  conductivity  of  electron  conductor 
particles  does  not  have  a  pronounce  effect  on  cell  performance, 
increasing  or  decreasing  the  ionic  conductivity  and  kinetic  con¬ 
stants  may  have  large  impact  on  polarization  resistance.  Best 
performance  has  been  found  to  fall  within  40-60  vol.%  of  elec¬ 
tron  conducting  particles. 

2.  Model  development 

Electrochemical  reaction  in  a  fuel  cell  composite  electrode 
occurs  at  the  contact  area  between  electronic  and  ionic  con¬ 
ducting  particles.  Furthermore,  the  electronic  conductor  (LSM 


Fig.  1.  Simplified  schematic  of  a  SOFC  cathode  showing  the  reaction  line  that 
starts  from  the  left  hand  side  of  the  electrode  and  extends  to  the  right  hand  side. 


for  the  cathode)  must  be  connected  contiguously  between  its 
contact  with  ionic  conductor  (typically  YSZ)  and  the  current 
collector  through  chains  of  other  electronic  conductors  to  be 
able  to  transfer  the  electron  to  the  reaction  point  and  similarly 
the  ionic  conductor  must  be  connected  to  the  electrolyte  to  con¬ 
duct  ions.  To  complete  the  process,  open  pores  must  be  present 
to  provide  reactants  for  the  reaction  and  to  carry  away  the  prod¬ 
ucts  of  the  electrochemical  reaction.  These  conditions  make  the 
active  site  or  active  three-phase  boundary  (TPB)  a  very  specific 
place  in  a  sense  that  not  all  three-phase  contact  regions  are  active 
as  assumed  by  the  current  models. 

An  ideal  structure  might  be  a  co-continuous  morphology  that 
is  commonly  used  in  the  polymer  composite  literature  wherein 
every  point  within  a  phase  is  topologically  connected  to  every 
other  point  within  the  same  phase.  Otherwise  islands  of  one 
phase  that  are  not  connected  to  others  are  inactive  and  do  not 
participate  in  the  electrochemical  processes.  Realization  of  such 
morphologies  depends  in  the  fabrication  processes  and  it  is  not 
commonly  found  in  SOFC  composite  electrodes. 

Fig.  1  illustrates  the  situation  for  a  simplified  2D  structure  of 
a  SOFC  cathode.  Black-coloured  particles  represent  LSM  par¬ 
ticles  while  YSZ  particles  are  represented  by  the  grey-coloured 
regions.  Supposing  that  pores  are  available  all  over  the  electrode 
to  transfer  reactant  (oxygen  in  case  of  cathode)  to  the  reaction 
sites,  it  is  clear  from  this  figure  that  not  all  contact  points  between 
LSM  and  YSZ  particles  are  electrochemically  active.  For  exam¬ 
ple  contact  point  4B’  is  not  electrochemically  active  because  the 
LSM  particle  is  not  connected  to  the  current  collector  while  the 
contact  point  ‘A’  is  active  as  a  result  of  capability  of  LSM  and 
YSZ  particles  to  conduct  electrons  and  ions,  respectively. 

Connecting  all  electrochemically  active  TPB’s  a  reactive  line 
can  be  constructed  for  this  2D  model  of  the  cathode  as  shown  in 
Fig.  1 .  The  length  of  this  line  which  is  a  measure  of  the  available 
sites  for  the  reaction  is  a  function  of  composition  of  electrode. 
It  is  believed  that  the  size  ratio  of  the  particles  also  has  great 
impact  on  the  length  of  this  line  [12,15]. 

Active  line  then  divides  the  electrode  into  two  parts,  the  upper 
part  which  is  basically  electron  conductor  and  the  lower  part 
which  is  ion  conductor.  All  YSZ  particles  falling  in  the  upper 
part  form  isolated  islands  which  are  neutral  to  electrons  and  LSM 
particles  below  the  active  line  are  neutral  to  ion  conduction  by 
the  same  token.  Although  the  isolated  particles  are  electrically 
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Fig.  2.  Simplified  schematic  of  a  SOFC  cathode  showing  the  reaction  line  that  starts  from  the  left  hand  side  of  the  electrode  and  extends  to  the  right  hand  side. 


neutral,  they  still  are  permeable  to  momentum,  mass  and  heat. 
Fig.  2  is  another  more  complex  geometry  obtained  from  ran¬ 
domly  mixing  of  equal  size  YSZ  and  LSM  particles  and  then 
transformation  of  the  geometry  into  a  continuous  media  of 
purely  electronic  or  ionic  conductors  with  pores  all  over  the  elec¬ 
trode.  Grey-coloured  continuous  domain  corresponds  to  YSZ 
while  white-coloured  continuous  media  indicates  LSM.  Note 
isolated  white-coloured  islands  all  over  the  electrode.  As  it  can 
be  seen  the  active  TPB  line  is  few  orders  of  magnitude  longer 
than  electrode  widths.  The  length  of  this  line  as  discussed  by 
the  other  people  [12,15]  is  dependent  on  the  ratio  of  electrode 
thickness  to  the  particle  size  and  also  depends  on  the  size  ratio 
of  electronic  and  ionic  conductor  particles. 

The  development  of  the  model  consists  of  the  following  steps: 

1 .  Creating  a  random  matrix  of  electronic  and  ionic  conductor 
particles  with  prescribed  composition  similar  to  Fig.  1. 

2.  Transforming  the  data  to  a  geometrical  form  similar  to  Fig.  2. 

3.  Model  the  transport  and  reaction  processes  in  each  domain 
using  any  multiphysics  solver. 

A  computer  code  has  been  written  to  carry  out  the  first  and 
second  steps  and  the  third  step  is  done  using  COMSOL®  a  com¬ 
mercial  finite  element  software. 

3.  Governing  equations 

Ohm’s  law  governs  the  electron  and  ion  transfer  inside  elec¬ 
tronic  and  ionic  conductors  along  with  the  charge  balance  equa¬ 
tion: 

V  (pel  =  pfjt  l  (1) 

V</>io  =  pjf/io  (2) 

/ei  =  —  Ao  =  i  (in  the  active  boundary)  (3) 

where  0  is  voltage,  yOeff  the  effective  resistivity  and  /is  the  current 
density.  Subscripts  ‘io’  and  ‘el’  represent  ionic  and  electronic 
conductors,  respectively.  Symbol  i  denotes  the  transfer  current 
density  of  the  cathode  which  is  an  indication  of  the  rate  of  elec¬ 
trochemical  reaction  taking  place  along  the  active  line  and  is 


expressed  by  the  classical  Butler- Volmer  equation: 


afri  _  e(l  -a)fr) 


where  io  is  the  exchange  current  density,  Co  the  oxygen  con¬ 
centration  and  Cq  is  the  equilibrium  oxygen  concentration. 
Overpotential  rj  in  this  equation  is  defined  as: 


rj  =  E  -  E o  =  (0io  -  0el)  -  E0 


where  Eo  is  equilibrium  potential  of  the  electrode  and  is  calcu¬ 
lated  from  well  known  Nemst  equation. 

Effective  porous  media  resistivity  for  electronic  and  ionic 
subdomain  is  calculated  from  pure  material  resistivity  by  the 
following  equation  [15]: 


Pk 

1  —  £ 


where  k  refers  to  the  electronic  or  ionic  conductor  and  s  is  the 
porosity  which  is  set  to  average  value  of  34.66%  [17]. 

The  diffusion  flux  through  the  cathode  can  be  determined 
using  the  well-known  Stefan-Maxwell  equation  [18]: 


n 


1 


-cVXi=  —^(xjNi-XiNj) 

j=hifj  7 


where  c  is  the  total  concentration  of  the  gas  mixture,  x  the  mole 
fraction  of  species  and  A/  is  the  diffusive  flux  of  specie  i.  is 
the  effective  binary  diffusion  coefficient  and  carries  the  effect  of 
both  ordinary  binary  diffusion  coefficient  and  Knudsen  diffusion 
coefficient  [19]: 


1 


D 


eff 

O2—  N2 


1  (  1 

n5,eff  nK,eff 

^o2-n2  uo2 


Binary  diffusion  coefficient  can  be  calculated  using  kinetic 
theory  by  Chapman-Enskog  relation  [18]: 


db  =  -mjmm 

lJ  Pcrn^D 


Knudsen  diffusion  coefficient  is  a  function  of  molecular 
velocity  and  radius  of  the  pores.  In  this  work  experimental  values 
reported  by  Zhao  et  al.  [17]  have  been  used. 


968 


A.  Abbaspour  et  al.  /  Journal  of  Power  Sources  161  (2006)  965-970 


Table  1 

Values  of  constants 


Parameter 

Value 

References 

Pel 

7.817  x  l0~5Qm 

[24] 

Pio 

0.442  Q  m 

[24] 

k) 

2000  A  m-2 

[19] 

T 

1073  K 

£ 

34.66% 

[17] 

X 

15.6 

[17] 

p.k,eff 

uo2 

0.624  cm2  s-1 

[17] 

Finally  to  calculate  effective  binary  diffusion  coefficients  for 
the  porous  media  from  the  values  obtained  using  Eq.  (8),  the 
following  equation  based  on  Cussler  [20]  has  been  used: 


r\B,eff  _  _  y^B 

Uij  ~  T  UH 


(10) 


Tortuosity,  r,  of  the  porous  media  is  set  to  15.6  [17].  Table  1 
summarizes  the  most  important  parameters  and  constants  used 
in  the  model. 


4.  Results  and  discussion 

Fig.  3  shows  the  results  obtained  from  the  model  for  a 
cathode  composed  of  about  2000  particles  in  a  mixture  of 
50vol.%  of  both  LSM  and  YSZ.  Applied  overpotential  to  the 
electrode  is  0.2  V.  It  can  be  seen  in  Fig.  3a  that  the  concen¬ 
tration  of  reactant  oxygen  is  not  uniform  along  the  cell  width. 
Instead  it  changes  depending  on  the  density  of  electrochemi- 
cally  active  points  and  the  minimum  concentration  occurs  near 
the  electrode-electrolyte  interface  at  regions  of  high  active  three- 
phase  boundary.  Fig.  3b  shows  the  variation  of  current  generation 
or  the  rate  of  electrochemical  reaction  along  active  boundary. 
This  figure  shows  that  although  the  concentration  of  oxygen  is 
lower  near  the  cathode-electrolyte  interface,  the  reaction  rate  is 
higher  in  this  region.  This  can  be  attributed  to  the  lower  ionic 
conductivity  of  YSZ.  Since  the  ionic  conductivity  of  YSZ  is 
much  smaller  than  the  electronic  conductivity  of  LSM  the  reac¬ 
tion  tends  to  occur  near  the  electrolyte  to  face  smaller  ionic 
resistance.  This  result  is  consistent  with  the  results  published  by 
Chan  and  Xia  [13]  who  reported  the  same  phenomena  for  the 
anode.  Although  this  trend  holds  for  the  practical  range  of  the 
value  of  parameters,  simulation  showed  that  one  cannot  always 
expect  higher  reaction  rates  near  the  electrode-electrolyte  inter¬ 
face.  In  fact  apart  from  activation  polarization  or  resistance, 
concentration  and  Ohmic  resistances  are  both  important  in  deter¬ 
mining  the  rate  of  electrochemical  reaction.  If  by  any  means 
oxygen  supply  to  the  reaction  sites  deep  inside  the  electrode 
decreases  it  would  result  in  electrode  starvation  at  that  points 
and  the  trend  will  be  altered.  This  can  happen  when  the  porosity 
of  electrode  is  not  enough  to  transport  the  oxygen  to  reaction 
sites  resulting  in  low  effective  diffusivity  and  thereby  starvation 
of  the  cell  near  the  electrode-electrolyte  interface.  Although 
it  is  not  the  case  for  LSM-YSZ,  high  electrocatalyst  activity 
resulting  in  high  exchange  current  density  can  also  result  in  the 
alteration  of  the  abovementioned  trend. 


Fig.  3.  (a)  Mole  fraction  of  oxygen  inside  the  SOFC  cathode  obtained 
by  the  model.  Upper  boundary  is  current  collector;  lower  boundary  is 
cathode-electrolyte  interface.  Mole  fraction  of  oxygen  in  the  current  collector 
side  is  set  to  0.21.  (b)  Current  generation  along  the  active  boundary. 


Fig.  4  shows  the  dimensionless  length  of  active  line  (length 
of  active  line  per  unit  length  of  electrode  width)  as  a  function 
of  LSM  volume  fraction.  As  pointed  out  by  most  of  the  authors 
[12,15]  the  length  of  active  line  or  the  active  area  assumes  its 


Fig.  4.  Effect  of  composition  on  the  electrode  active  area. 
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LSM  volume  fraction 

Fig.  5.  Effect  of  composition  on  the  overall  polarization  resistance  of  the  elec¬ 
trode. 

maximum  value  when  the  electrode  is  made  up  of  equal  amounts 
of  electronic  and  ionic  conductors  of  the  same  size. 

The  effect  of  electrode  composition  on  the  overall  electrode 
resistance  is  shown  in  Fig.  5.  Overall  polarization  resistance  is 
defined  as: 

ii 

R=j  (11) 

where  I  is  the  total  current  density  of  the  electrode.  As  we  can 
see  the  polarization  resistance  for  the  cell  is  minimum  some¬ 
where  around  50-60  vol.%  of  the  electron  conductor.  Although 
the  length  of  active  line  is  maximum  at  exactly  50%  and  the  vari¬ 
ations  are  symmetric  according  to  Fig.  4,  for  polarization  it  is  not 
symmetric  and  minimum  resistance  does  not  happen  in  the  same 
composition.  Even  though  activation  polarization  is  minimum 
when  the  length  of  active  line  is  maximum  (50%),  total  polar¬ 
ization  resistance  depends  on  some  other  factors  which  tend  to 
affect  it.  While  the  length  of  active  area  is  only  a  function  of  the 
geometry  of  the  electrode,  the  polarization  resistance  depends 
not  only  on  the  geometry  but  also  on  other  parameters  such  as 
conductivities,  diffusivity  and  electrode  porosity.  That  is  why 
different  researchers  obtained  slightly  different  optimum  com¬ 
position  of  the  electrode  depending  on  the  operating  condition 
and  parameters  of  their  simulations  or  experiments  while  every¬ 
body  agrees  that  the  maximum  active  line  or  area  occurs  when 
equal  amounts  of  LSM  and  YSZ  particles  of  the  same  size  are 
used. 

Kenjo  and  Nishiya  [3]  and  Haanappel  et  al.  [21]  reported  the 
best  performance  for  the  electrode  for  the  weight  ratio  of  LSM 
and  YSZ  equal  to  1.  While  the  first  group  did  not  mention  the 
size  of  the  particles,  second  group  used  1  p,m  LSM  particles.  On 
the  other  hand,  Ostergard  et  al.  [7]  and  Kim  et  al.  [11]  obtained 
the  minimum  polarization  resistance  in  their  experiments  when 
YSZ  content  was  around  40%.  Juhl  et  al.  [22]  reported  the  same 
results.  On  the  simulation  side  also  there  are  differences  between 
reported  data.  While  Costamagna  et  al.  [12]  report  the  minimum 
polarization  resistance  at  about  37%  of  electronic  conductor, 
Sunde  [23]  and  Chan  et  al.  [14]  believe  that  for  a  wide  range  of 


Fig.  6.  Parametric  study  of  the  effect  of  diffusion  coefficient  on  the  overall 
polarization  resistance  of  the  electrode. 

electrode  composition  of  about  40-60%  the  polarization  resis¬ 
tance  remains  constant  at  its  minimum  value. 

Figs.  6  and  7  show  a  parametric  study  of  the  effect  of  dif¬ 
fusion  coefficient  and  the  conductivity  of  ion  conductor  on  the 
performance  of  the  electrode.  Decreasing  the  diffusion  coeffi¬ 
cient  which  could  be  a  result  of  decreasing  porosity,  can  result  in 
diffusion  controlled  process.  According  to  the  model  increasing 
the  volume  fraction  of  electronic  conductor  pushes  the  elec¬ 
trochemical  reaction  sites  towards  the  electrolyte,  far  from  the 
air  channel.  This  in  turn  increases  the  concentration  polariza¬ 
tion  resistance  and  for  small  values  of  diffusion  coefficient  as 
it  is  shown  in  Fig.  6  concentration  polarization  resistance  out¬ 
weighs  the  increase  in  the  length  of  active  line  or  decrease  in 
activation  polarization  resistance.  Fig.  7  is  a  parametric  study 
on  the  effect  of  conductivity  of  ion  conductor  in  the  cathode. 
Since  the  conductivity  of  ion  conducting  media  is  usually  very 
low,  change  in  the  conductivity  has  pronounced  effect  on  the 
performance  of  the  electrode.  This  change  in  conductivity  can 
be  applied  through  changing  the  geometric  features  like  parti¬ 
cle  size  and  porosity  or  variations  in  operating  conditions  like 
temperature.  For  high  values  of  conductivity  the  minimum  of 
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LSM  volume  fraction 


Fig.  7.  Parametric  study  of  the  effect  of  the  conductivity  of  ionic  conductor  on 
the  overall  polarization  resistance  of  the  electrode. 
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polarization  resistance  shifts  exactly  to  50%vol.  of  each  com¬ 
ponent  and  the  shape  of  the  curve  becomes  more  symmetric 
just  like  active  area  or  active  length  curve  of  Fig.  4.  The  reason 
is  that  increasing  the  conductivity  decreases  the  Ohmic  resis¬ 
tance  and  thereby  activation  polarization  determines  the  whole 
trend  of  the  polarization  resistance  curve.  Activation  polariza¬ 
tion  is  minimum  when  the  active  area  is  maximum  or  when  the 
composition  is  50%.  A  slight  increase  in  polarization  resistance 
with  increasing  the  LSM  volume  fraction  in  the  first  stages  is 
due  to  the  fact  that  by  increasing  LSM  composition  in  the  first 
stages  we  are  just  introducing  some  small  isolated  patches  of 
LSM  into  YSZ  which  decreases  the  conductivity  of  YSZ  media 
hence  increasing  the  polarization  resistance.  Same  observation 
is  reported  for  conductivity  by  Sunde  [23]. 

5.  Conclusions 

A  new  2D  model  based  on  random  packing  of  electronic 
and  ionic  particles  has  been  proposed.  The  model  of  transport 
processes  based  on  this  geometrical  approach  shows  very  good 
agreement  with  the  other  observed  experimental  and  simulation 
data.  This  framework  offers  its  own  values  and  advantages,  since 
in  this  model  we  can  identify  every  reaction  point,  and  hence  it 
gives  a  better  image  of  composite  electrode  performance.  It  can 
be  used  to  study  the  reaction  kinetics  of  electrode  reaction. 

We  must  extend  this  to  3D  framework,  which  remains  a  purely 
geometrical  problem  to  be  solved,  although  the  extra  degree 
of  freedom  offered  in  a  3D  geometry  can  alter  the  topological 
structure  of  the  TPB  significantly.  Using  this  framework,  we  can 
also  investigate  the  effectiveness  of  several  alternate  structured 
or  patterned  topologies  to  assess  the  most  promising  patterns. 
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